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Abstract 



We numerically study in the one-dimensional case the validity of the func- 
tional calculated by Graham and coworkers (R. Graham and T. Tel, Phys. 
Rev. A 42, 4661 (1990), O. Descalzi and R. Graham, Z. Phys. B 93, 509 
(1994)) as a Lyapunov potential for the Complex Ginzburg-Landau equation. 
In non-chaotic regions of parameter space the functional decreases monoton- 
ically in time towards the plane wave attractors, as expected for a Lyapunov 
functional, provided that no phase singularities are encountered. In the phase 
turbulence region the potential relaxes towards a value characteristic of the 
phase turbulent attractor, and the dynamics there approximately preserves a 
constant value. There are however very small but systematic deviations from 
the theoretical predictions, that increase when going deeper in the phase tur- 
bulence region. In more disordered chaotic regimes characterized by the pres- 
ence of phase singularities the functional is ill-defined and then not a correct 
Lyapunov potential. 
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I. INTRODUCTION 



The Complex Ginzburg-Landau Equation (CGLE) is the amplitude equation describing 
universal features of the dynamics of extended systems near a Hopf bifurcation Jl]PJ] . 

d t A = aA + (D r + iDi)V 2 A - (b r + ih) \ A | 2 A . (1.1) 

Examples of this situation include binary fluid convection ||, transversally extended lasers 
0] and chemical turbulence [[J. We will considered here only the one- dimensional case, 
A = A(x,t), with x G [0,L]. Suitable scaling of the complex amplitude A, space, and time 
shows that for fixed sign of a there are only three independent parameters in ( |1 . 1|) (with 
D r and b r > that we assume henceforth). They can be chosen to be L, c\ = Di/D r , and 
c 2 = bi/b r . 

The CGLE for a > displays a rich variety of complex spatio-temporal dynamical 
regimes that have been recently classified in a phase diagram in the parameter space {ci, C2} 
|6|-[S[]. It is commonly stated that such nontrivial dynamical behavior, occurring also in other 
nonequilibrium systems, originates from the non-potential or non-variational character of 
the dynamics [[5J. This general statement needs to be qualified because it involves some 
confusion in the terminology. For example the term "non-variational" is often used meaning 
that there is no Lyapunov functional for the dynamics. But Graham and co-workers, in a 
series of papers [|I0|-|T4||, have shown that a Lyapunov functional does exist for the CGLE, 
and they have constructed it approximately in a small-gradient approximation. The correct 
statement for the CGLE is that it is not a gradient flow. This means that there is no real 
functional of A from which the right hand side of ( |1 . 1|) could be obtained by functional 
derivation. 

Part of the confusion associated with the qualification of "nonvariational" dynamics 
comes from the idea that the dynamics of systems having non-trivial attractors, such as 
limit cycles or strange chaotic attractors, can not be deduced from the minimization of a 
potential which plays the role of the free energy of equilibrium systems. However, such 
idea does not preclude the existence of a Lyapunov functional for the dynamics. The Lya- 
punov functional can have local minima which identify the attractors. Once the system has 
reached an attractor which is not a fixed point, dynamics can proceed on the attractor due 
to "nonvariational" contributions to the dynamical flow which do not change the value of the 
Lyapunov functional. This just means that the dynamical flow is not entirely determined 
once the Lyapunov functional is known. This situation is very common and well known in 
the study of dynamical properties within the framework of conventional statistical mechan- 
ics: The equilibrium free energy of the system is a Lyapunov functional for the dynamics, 
but equilibrium critical dynamics |15| usually involves contributions, such as mode-mode 



coupling terms, which are not determined just by the free energy. The fact that the dynam- 
ical evolution is not simply given by the minimization of the free energy is also true when 
studying the nonequilibrium dynamics of a phase transition in which the system evolves 
between an initial and a final equilibrium state after, for example, a jump in temperature 
across the critical point [ |T6|j . 

A Lyapunov functional plays the role of a potential which is useful in characterizing 
global properties of the dynamics, such as attractors, relative or nonlinear stability of these 
attractors, etc. In fact, finding such potentials is one of the long-sought goals of nonequilib- 
rium physics [|P71 , |l8f , the hope being that they should be instrumental in the characterization 
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of nonequilibrium phenomena through phase transitions analogies. The use of powerful and 
very general methods based on these analogies has been advocated by a number of authors 
T9| , |20| j6|-p| . In this context, it is a little surprising that the finding of a Lyapunov functional 



for the CGLE [|12l-(14}| has not received much attention in the literature. A possible reason 
for this is that the construction of nonequilibrium potentials has been historically associated 
with the study of stochastic processes, in particular in the search of stationary probability 
distributions for systems driven by random noise |T^,|TB|,^T| . We want to make clear that 



the finding of the Lyapunov functional for the CGLE fl!q-|14|, as well as the whole approach 
and discussion the present paper is completely within a purely deterministic framework and 
it does not rely on any noise considerations. A second possible reason for the relative little 
attention paid to the Lyapunov functional for the CGLE is the lack of any numerical check 
of the uncontrolled approximations made on its derivation. The main purpose of this paper 
is precisely to report such numerical check of the results of Graham and collaborators, thus 
delimiting the range of validity of the approximations involved. We also provide a charac- 
terization of the time evolution of the Lyapunov functional in different regions of the phase 
diagram of the CGLE which illustrates the use of such potential. 

Our main findings are that the expressions by Graham and coworkers behave to a good 
approximation as a proper Lyapunov potential when phase singularities (vanishing of the 
modulus of A) are not present. This includes non-chaotic regimes as well as states of phase 
turbulence. In this last case some small but systematic discrepancies with the predictions 
are found. In the presence of phase singularities the potential is ill-defined and then it is 
not a correct Lyapunov functional. 

The paper is organized as follows. For pedagogical purposes, we first discuss in Sect. II 
a classification of dynamical flows in which notions like relaxational or potential flows are 
considered. The idea of a potential for the CGLE is clearer in this context. In Sect. Ill we 
review basic phenomenology of the CGLE and the main analytical results for the Lyapunov 
functional of the CGLE. Sections IV and V contain our numerical analyses. Section IV is 
devoted to the Benjamin- Feir stable regime of the CGLE and Sect. V to the Phase Turbulent 
regime. Our main conclusions are summarized in Sect. VI. 



II. A CLASSIFICATION OF DYNAMICAL FLOWS 

In the following we review a classification of dynamical systems that, although rather well 
established in other contexts [ |TT| , |r8| , it is often overlooked in general discussions of deter- 
ministic spatio-temporal dynamics. Non-potential dynamical systems are often defined as 
those for which there is no Lyapunov potential. Unfortunately, this definition is also applied 
to cases in which there is no known Lyapunov potential. To be more precise, let us consider 
dynamical systems of the general form 

d t A t = Vi[A\ (2.1) 

where A{ represents a set of, generally complex, dynamical variables which are spatially 
dependent fields: Ai = Aj(x, t). Vi[A] is a functional of them. The notation A* represents 
the complex conjugate of Aj and for simplicity we will keep the index i implicit. Let us now 
split V into two contributions: 
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V[A] = G[A] + N[A] , 



(2.2) 



where G, the relaxational part, will have the form 

with F a real and scalar functional of A. V is an arbitrary hermitic and positive-definite 
operator (possibly depending on A). In the particular case of real variables there is no need 
of taking the complex conjugate, and hermitic operators reduce to symmetric ones. The 
functional N[A] in ( J2.2| ) is the remaining part of V[A\. The important point is that, if 
the splitting ( |2.2| ) can be done in such a way that the following orthogonality condition is 
satisfied (c.c. denotes the complex conjugate expression): 

rfx (^|^[A(x)] + c.c.) = , (2.4) 

then the terms in N neither increase nor decrease the value of F, which due to the terms in 
G becomes a decreasing function of time: 

dF W*'"' < . (2.5) 
at 

If F is bounded from below then it is a Lyapunov potential for the dynamics ( |2.1| ). Equation 
(UD with N = V - G, that is 



can be interpreted as an equation for the Lyapunov potential F associated to a given dynam- 
ical system (|2.1| ). It has a Hamilton- Jacobi structure. When dealing with systems perturbed 
by random noise, V is fixed by statistical requirements, but in deterministic contexts such 
as the present paper, it can be arbitrarily chosen in order to simplify ( [2.6|) . 

Solving ( |2.6| ) is in general a difficult task, but a number of non-trivial examples of the 
splitting (|2.2|) - (|2.5|) exist in the literature. Some of these examples correspond to solutions 
of ( |2.6| ) found in the search of potentials for dynamical systems [fT^ , P~D|JTT|1 . Other examples 
just correspond to a natural splitting of dissipative and non-dissipative contributions in 
the dynamics of systems with well established equilibrium thermodynamics, as for example 
models of critical dynamics |15| or the equations of nematodynamics in liquid crystals ||22| . 



Once the notation above has been set-up, we can call relaxational systems those for 



which there is a solution F of (2J3) such that N = 0, that is all the terms in V contribute to 
decrease F. Potential systems can be defined as those for which there is a nontrivial (i.e. a 
non-constant) solution F to ( |2.6| ). In relaxational systems there is no long-time dynamics, 
since there is no time evolution of A once a minimum of F is reached. On the contrary, for 
potential systems for which N ^ 0, the minima of F define the attractors of the dynamical 
flow, but once one of these attractors is reached, nontrivial sustained dynamics might exist 
on the attractor. Such dynamics is determined by iV and maintains a constant value for the 
functional F. 

A possible more detailed classification of the dynamical flows is the following: 
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1.- Relaxational gradient flows: Those dynamical systems for which N = with T propor- 
tional to the identity operator. In this case the time evolution of the system follows the 
lines of steepest descent of F. A well known example is the so called Fisher-Kolmogorov 
equation, also known as model A of critical dynamics , or (real) Ginzburg-Landau 
equation for a real field A(x, t): 



A = aA + -iV 2 A- (3 \ A | 2 A 



(2.7) 



where a, 7 and (5 are real coefficients. This equation is of the form of Eqs. ( ^.1| )-( P75D 
with N — 0, r = 1, and F = F GL [A], the Ginzburg-Landau free energy: 



F GL [A] 



dx —a I A 1 2 +7 



VA \ 2 +-\A 



(2.8) 



2.- Relaxational non-gradient flows: Still N = but with T not proportional to the iden- 
tity, so that the relaxation to the minimum of F does not follow the lines of steepest 
descent of F. The matrix operator V might depend on A or involve spatial deriva- 
tives. A well known example of this type is the Cahn-Hilliard equation of spinodal 
decomposition, or model B of critical dynamics for a real variable A. [ D| : 



The symmetric and positive-definite operator (—V 2 ) has its origin in a conservation 
law for A. 

3.- Non-relaxational potential flows: N does not vanish, but the potential F, solution 
of ( [2.6D exists and is non-trivial. Most models used in equilibrium critical dynamics 
[y| include non-relaxational contributions, and therefore belong to this category. A 



particularly simple example is 

^ = -(1 + 0%^, (2.10) 

where now A is a complex field. Notice that we can not interpret this equation as 
being of type 1, because (1 +i) is not a hermitic operator, but still Fql is a Lyapunov 
functional for the dynamics. Equation (|2.10|) is a special case of the Complex Ginzburg- 
Landau Equation (CGLE), in which V[A] is the sum of a relaxational gradient flow 
and a nonlinear-Schrodinger-type term N[A] = —i SF f^l A ^ ■ 

The general CGLE is of the form (|2.7|) but A is complex and a, 7 and (3 are 
arbitrary complex numbers. For the special case in which j^A = > as f° r example 
in (|2.10|) , the Lyapunov functional for the CGLE is known exactly p3| . Such choice of 
parameters has important dynamical consequences ||24j|. Beyond such special cases, the 



calculations by Graham and coworkers indicate p3] , p^| that the CGLE, a paradigm 
of complex spatio-temporal dynamics, might be classified within this class of non- 
relaxational potential flows because a solution of ( |2.6| ) is found. The difficulty is that 



the explicit form of the potential is, so far, only known as a uncontrolled small-gradient 
expansion. 
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4.- 



Non-potential flows: Those for which the only solutions F of (|2.6|) are the trivial ones 
(that is F = constant). Hamiltonian systems as for example the nonlinear Schrodinger 
equation are of this type. 



III. A LYAPUNOV FUNCTIONAL FOR THE CGLE 



It is well known that for a < the one dimensional CGLE (O) has A = as a stable 



solution, whereas for a > there are Travelling Wave (TW) solutions of the form 

A k = A e i{kx+UJt)+V0 (3.1) 



with A = y (a — D r k 2 )/b r , \k\ < Ja/b r , and u = (6^a + D_k 2 )/b r . We have introduced 

D- = DA - DA . (3.2) 

(fo is any arbitrary constant phase. 

The linear stability of the homegeneus solution ( (|3.1| ) with k = ) with respect to long 
wavelength fluctuations divides the parameter space {c\,C2} in two regions: the Benjamin- 



Feir (BF) stable and the BF unstable zone. This line is given by p5jp6 



D+ = D r b r + DA = 0, (3.3) 

In the BF unstable region (D + < 0) there are no stable TW solutions, while in the BF 
stable region (D + > 0) TWs with a wavenumber k < k^ are linearly stable. For k > 
kg, TWs become unstable through the long wavelength instability known as the Eckhaus 
instability [^3,^]. The Eckhaus wavenumber Ue is given by 

2 ab r D + 

E D r {3D + b r + 2D_bi) 1 1 



Recent numerical work for a > and L large |Q-^,29| has identified regions of the 



parameter space displaying different kinds of regular and spatio-temporal chaotic behavior 
(obtained at long times from random initial conditions and periodic boundary conditions), 
leading to a "phase diagram" for the CGLE. The five different regions, each leading to a 
different asymptotic phase, are shown in Fig. |I] as a function of the parameters c\ and C2 
(a > 0, L large). Two of these regions are in the BF stable zone and the other three in the 
BF unstable one. One of the main distinctions between the diferent asymptotic phases is in 
the behavior of the modulus of A at long times. In some regions it never vanishes, whereas 
in others it vanishes from time to time at different points. A more detailed description of 
the asymptotic behavior in the different regions is as follows: 

1. Non-Chaotic region. The evolution here ends in one of the Eckhaus-stable TW solu- 
tions for almost all the initial conditions. 

2. Spatio- Temporal Intermittency region. Despite the fact that there exist stable TW, 
the evolution from random initial conditions is not attracted by them but by a chaotic 
attractor in which typical configurations of the field A consist of patches of TW inter- 
rupted by turbulent bursts. The modulus of A in such bursts typically touches zero 
quite often. 



6 



3. Defect Turbulence. This is a strongly disordered phase in which the modulus of A 
has a finite density of space-time zeros. In addition the space and time correlation 
functions have a quasi-exponential decay |||7[]. 

4. Phase Turbulence. This is a weakly disordered phase in which \A(x, t)\ remains away 
from zero. The temporal correlations decay slower than exponentially |].[7| . 

5. Bi-Chaos region. Depending on the particular initial condition, the system ends on 
attractors similar to the ones in regions 3, 4, or in a new attractor in which the 
configurations of A consists of patches of phase and defect turbulence. 

An approximate Lyapunov functional for the CGLE was calculated by Graham and 
collaborators |l~3|,^4|,p0| . Earlier attempts to find a Lyapunov functional were based on 



polynomial expansions p3| , |31| -p^1 , while more recent and successful approaches focussed in 
solving the Hamilton- Jacobi equation (|2.6|) with T = 1 in different ways. This was done 
first by a minimization procedure involving an action integral JTTJ-12], and more recently by 
a more direct expansion method |T3|JI4]j30[1 . This last method provides also expressions in 
higher dimensions, but we will restrict here to the one- dimensional case. In any case, the 
solution involves an uncontrolled gradient expansion around space-independent solutions of 
the CGLE. Such expansion obviously limits the validity of the result to regions in the phase 
diagram in which there are not strong gradients. Since the expansion was actually performed 
in polar coordinates, this excludes the regions in which zeros in the modulus of A are typical, 
since the phase of A becomes singular there. In particular Spatio-temporal intermittency 
regimes, Bi-chaos and Defect Turbulence are out of the range of validity of Graham's expan- 
sion. The meaningfulness of the potential in the other regions of parameter space remains 
still an open question because of the uncontrolled small gradient approximations used to 
calculate it, and calls for some numerical check. 

In their solution of the Hamilton Jacobi equation, Graham and collaborators find dif- 
ferent branches of the Lyapunov functional with expressions valid for different values of the 
parameters. In particular they identify the BF line ( |3.3|) as separating two branches of the 
solution to (|2.6|). 

The explicit expressions (obtained with T = 1) are given in polar coordinates: 



A(x, t) = r(x, t)e 



iip(x,t) 



(3.5) 



In terms of the amplitude r, the phase <p, and their spatial derivates (denoted as r x , (p x ,ip xx , 
etc.) the Lyapunov functional per unit of length $ = F/L was found OJTl], for a < 0: 



b r r A - 2ar 2 + 2 



D r + 



D.bir 4 
3(a — b r r 2 ) 2 



2D_r 3 

r " ~ * u 2S *<P* + 2D rr 2 <fil \dx (3.6) 
Ma — b r r z ) 



We note that even in this relatively simple case a < 0, the result for $ is only approximate 
and its structure reveals a highly non-trivial dynamics. 

For a > 0, in the BF stable region (D + > 0) the expression for $ results: 



by 



2ar' 



+ 



(A ir + Bjryi + {A 2 r + B 2 /r)r x <p x + 2{D T r 2 - D_b t a/ | b \ 2 b r )<?l 
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D_D r bi 4 / D 2 _a D— . , nr0l A 2 2D_D r 2 2 



+ _ 2bl) ^ - ADAD ~ (l + Hby/a) )r x ^ 

3bf | b | 2 1 r r 2 36^?r V 1 — b r r 2 /aj " ' " ' ' 



jrfx (3.7) 



where 



A 1 = 2(D r + biD_/3b 
A 2 = 2D_/b 



V ) ■ 



2D-6<o , n 2 2 



(3.8) 



Clearly, $ is ill-defined when r = 0. 

By writing-out the Euler-Lagrange equations associated to the minimization of $ the TW 
solutions ( |3.1|) are identified as local extrema of $. Since they occur in families parametrized 
by the arbitrary phase ipo, the minima associated to the TW of a given k are not isolated 
points but lay on a one-dimensional closed manifold. The non-variational part of the dy- 
namics (N in Q2.2D ) can be explicitly written-down by substracting G = —\4rz with F = L$ 
to the right-hand-side of (|1.1|) . It is seen to produce, when evaluated on the manifold of 
minima of $ with a given k, constant motion along it. This produces the periodic time 
dependence in ( |3.1| ) and identify the TW attractors as limit cycles. 

The value of k for which the corresponding extrema change character from local minima 
to saddle points is precisely the Eckhaus wavenumber fog. It is remarkable that, although 
expression ( |3.7| ) was obtained in a gradient expansion around the homogeneous TW, their 
minima identify exactly all the TWs of equation ( p. . 1|) , and their frequencies and points of 
instability are also exactly reproduced. This gives confidence on the validity of Graham's 
approximations. It should be stressed however that they are not exact and can lead to 
unphysical consequences. For instance, the value of the potential $ evaluated on a TW of 



wavenumber k (\k\ < Ja/b r ) is JT2 



^ = $[^] = ^ 2 (l-^)+^ (3.9) 

where §k=o — —a 2 /b r . For a range of parameter values this expression gives mathematical 
sense to the intuitive fact that the closer to zero is k the more stable is the associated TW 
(because its potential is lower). But for some parameter values the minimal potential cor- 



responds to large wavenumbers close to ±ya/b r . This is counterintuitive and calls for some 
numerical t est. T he test will be described below and it will be shown that the wavenumbers 
close to ±^Ja/b r are out of the range of validity of the small gradient approximations leading 

to©. 

We already mentioned in the previous section that the Lyapunov functional for the 
CGLE is exactly known for special values of the parameters [T2"|,[T3"|.|2~3| . This happens for 
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D- = D r b{ — Dib r = 0, which lies in the BF-stable region as indicated in Fig. |l] . In this 
case it is clear that (11 . II) can be written as 



.4 



I SF GL [A] 
2 5 A* 



+ ibi 



A | 2 +^V 2 

b r 



A 



(3.10) 



where .Fg£,L4] is (|2.S| ) for complex A and with a = 2a, f3 = 2b r , and 7 = 2D r . It is readily 
shown that the term proportional to 6j is orthogonal to the gradient part, so that Fql is an 
exact solution of (|2.6|) for these values of the parameters, and ( 3.10|) is a relaxational non- 
gradient flow (see classification in section ||). It is seen that the approximate expressions 

leading both to the same expression: 

(3.11) 



1)D and (|3.7| ) greatly simplify when D- 

L$ = J |-2ar 2 + b r r 4 + 2D T rr\ + 2D r r 2 <fl\dx 



When expressed in terms of A and A* it reproduces Fql in ( |3.10| ). Thus the gradient 
expansion turns out to be exact on the line _D_ = 0. 

In the Benjamin-Feir unstable region (a > 0, D + < 0) the gradient expansion for $ 
becomes |Oj30|: 



+ 



b r r 4 — 2ar 2 + 

b r 



{A x r + #i/r 2 )r 2 + (A 2 r + B 2 /r)r x <f x + 2D r {r 2 



br 



2a?r 2 



B 2 AD 2 a 2 



4 



l> 2 . 



a 



A 2/ A 2 r 2 b r , 
+ Di) ln( 



1 + 



D r \b\ 2 +2b r D. 



biDjl 



b r r 2 



ln( 



,r 2 b r , 



2b r 

^Px^x^xx 



4 



+ 



2 

Vxx 



where, in addition to the previous definitions 

2D-bi,a 



Dib r )v 2 x r 2 x 



dx 



(3.12) 



B 2 



m 1 b 1 2 



n-b% 



~(D r bi + D { b r 
bi 



(3.13) 

It was noted before that this expression can be adequate, at most, for the Phase Turbulent 
regime, since in the other BF unstable regimes \A\ vanishes at some points and instants, so 
that ( pT2D is ill-defined. 

The long time dynamics occurs in the attractor defined by the minima of $. The Euler- 
Lagrange equations associated to the minimization of ( |3.12| ) lead to a relationship between 
amplitude and phase of A which implies the well known adiabatic following of the amplitude 
to the phase dynamics commonly used to describe the phase turbulence regime by a nonlinear 
phase equation. The explicit form of this relationship is 



a D r ^ 2 Di 2 biD 2 4 , n D r DA 



br b r 1 y) b r * 2abl 



vv + 2- 



abl 



b-D 2 

2^Vy>V W + 
abi 



DrDjbi 

abl 



D 



ab r 



(vV) s 



(3.14) 
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It defines the attractor characterizing the phase turbulent regime. Dynamics in this 
attractor follows from the nonrelaxational part N in 



r2|) . When ( 3.14 ) is imposed in 
such nonrelaxational part of the dynamics the generalized Kuramoto-Shivashinsky equation 
containing terms up to fourth order in the gradients |34| is obtained fl4] , |30"|l . 



We finally note that in the phase turbulent regime the Lyapunov functional $ gives the 
same value [^IjpOf when evaluated for any configuration satisfying (|3.14|) , at least within 
the small gradient approximation. This corresponds to the evolution on a chaotic attractor 
(associated to the Kuramoto-Sivashinsky dynamics coming from N) which is itself embedded 
in a region of constant <3? (the potential plateau [IS]). This plateau consists of the functional 
minima of $ ( |3.14| ). All the (unstable) TW are also contained in the same plateau, since 
they satisfy fTUl) . 



IV. NUMERICAL STUDIES OF THE LYAPUNOV FUNCTIONAL IN THE 
BENJAMIN-FEIR STABLE REGIME 

We numerically investigate the validity of <&[A] in ( |3.6j ), fl3.7| ), and ( |3.12| ) as an approx- 
imate Lyapunov functional for the CGLE. When evaluated on solutions A(x, t) of ( |1 . 1|) 
it should behave as a monotonously decreasing function of time, until A(x, t) reaches the 
asymptotic attractor. After then, $ should maintain in time a constant value characteristic 
of the particular attractor. 

All the results reported here were obtained using a pseudo-spectral code with periodic 
boundary conditions and second-order accuracy in time. Spatial resolution was typically 
512 modes, with runs of up to 4096 modes to confirm the results. Time step was typically 
At = .1 except when differently stated in the figure captions. Since very small effects have 
been explored, care has been taken of confirming the invariance of the results with decreasing 
time step and increasing number of modes. System size was always taken as L = 512, and 
always D r = 1 and bi = —1, so that c\ = Di and C2 = —l/b r . When a random noise of 
amplitude e is said to be used as or added to an initial condition it means that a set of 
uncorrelated Gaussian numbers of zero mean and variance e 2 was generated, one for each 
collocation point in the numerical lattice. 



A. Negative a 

The uniform state A = is stable for a < 0. We start our numerical simulation with 
a plane wave A = Aoe tkx of arbitrary wavenumber k = 0.295 and arbitrary amplitude 
Aq = 1 (note that the TWs fl3.1|) do not exist for a < 0), and calculate $ for the evolving 
configurations. In order to have relevant nonlinear effects during the relaxation towards 
A = we have chosen a small value for the coefficient of the linear term (a = —0.01). The 
remaining parameters were Di = 1 and b r = 1.25 (ci = 1, c% — —0.8). Despite the presence 
of non-relaxational terms in (|l.lf) , $ decreases monotonously (see Fig. [|) to the final value 
$(t = oo) = $L4 = 0] = confirming its adequacy as a Lyapunov potential. 
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B. Positive a. Benjamin-Feir stable regime 



We take in this section always a — 1. Non-chaotic (TW) states and Spatio-Temporal 
Intermittency are the two phases found below the BF line in Fig. [I]. We first perform several 
numerical experiments in the non-chaotic region: 

A first important case is the one on the line _D_ = 0, for which ( |3.11| ) is an exact 



Lyapunov functional F GL . We take Di = — 1 and b r = 1 (ci = c 2 = — 1), on the D_ = 
line, and compute the evolution of $ = along a solution of ( |1 . 1| ) , taking as initial 
condition for A a Gaussian noise of amplitude e = 0.01. Despite of the strong phase gradients 
present specially in the initial stages of the evolution, and of the presence of non-relaxational 



terms, $ decays monotonously in time (Fig. |VI1| ). The system evolved towards a TW 
attractor of wavenumber k = 0.0245. The value of $ in such state is, from Eq. ( |3.9| ), 
$fc =0 .0245 = —0.998796. It is important to notice that our numerical solution for A and 
numerical evaluation of the derivatives in <3> reproduce this value within a 0.3% in the 
last time showed in Fig. reffig3, and continues to approach the theoretical value for the 
asymptotic attractor at longer timesQ. 

We continue testing the Lyapunov functional for = 1, b r = 1.25 (ci = 1 ,C2 = —0.8. 
This is still in the non-chaotic region but, since _D_ ^ 0, $ is not expected to be exact, but 
only a small gradient approximation. We check now the relaxation back to an stable state 
after a small perturbation. As initial condition we slightly perturb a TW of Eckhaus-stable 
wavenumber (k = 0.13 < /eg) by adding random noise of amplitude e = 0.09. <£> decays 
monotonously (Fig. |[) from its perturbed value to the value <£> fc=0 .i3 = —0.796632 as the 
perturbation is being washed out, as expected for a good Lyapunov functional. 

A more demanding situation was investigated for Di = —1 and b r = 0.5 (again in the 
non-chaotic region, c\ = —1 and C2 = —2, and D^ ^ 0). Two TW of different wavenumbers 
(ki = 0.4, £>2 = 0.08, both Eckhaus-stable) were joined and the resulting state (see inset in 
Fig. ^|) was used as initial condition. The TW of smaller wavenumber advances into the 
other, in agreement with the idea that it is nonlinearly more stable since it gives a smaller 
value to the potential. As the difference between the two frequencies is large the speed at 
which one wave advances onto the other is quite large. The interface between the two TWs 
contains initially a discontinuity in the gradient of the phase which is washed out in a few 
integration steps. An important observation is that during the whole process the modulus 
of A(x, t) never vanishes and then the winding number, defined as 



L 







Vipdx (4.1) 



remains constant [y = 20) (with periodic boundary conditions v is constant except at the 
instants in which the phase becomes singular, that is when r = 0). After the TW with the 
smallest wavenumber completely replaced the other, still a phase diffusion process in which 
the wave adjusts its local wavenumber to the global winding number occurs. The state (limit 



If a smaller time step is used greater accuracy is obtained. For example, if the time step is 
reduced to 0.05 the value of $ is reproduced within 10 _7 %. But this takes quite a long computing 
time. 
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cycle) finally reached is a TW of k = 2ttv/ L = 0.245. Despite of the complicated and non- 
relaxational processes occurring $ behaves as a good Lyapunov functional monotonously 
decreasing from the value $(£ = 0) = —1.825 corresponding to the two-wave configuration 
to the value $ = —1.863 of the final attractor (Fig. reffig5). It would be interesting, as 



happens in some relaxational models ||35|| , finding some relationship between the speed of 
propagation of the more stable wave onto the less stable one and the difference in $ between 
the two states. 

The good behavior of $ will be obviously lost if the field A(x, t) vanishes somewhere 
during the evolution. As the next numerical experiment (for Di = 1 and b T = 1.25, that 
is C\ = 1, Ci = —0.8) we used as initial condition a small (e = 0.01) random Gaussian 
noise. The system was left to evolve towards its asymptotic state (a TW). Fig. |] shows 
that after a transient $ monotonously decreases. During the initial transient it widely 
fluctuates, increasing and decreasing and loosing then its validity as a Lyapunov functional. 
This incorrect behavior occurs because during the initial stages A(x,t) is small and often 
vanishes, changing v. When A (and then r) vanishes the phase and ( |3.7|) are ill-defined 
and out of the range of validity of a small gradient approximation. Note the contrast with 
the case D- = in which the potential is exact and well behaved even when v is strongly 
changing. The particular values of the maxima and minima during the transient in which 
v is changing depend on the spatial and temporal discretization, since it is clear from fl3.7|) 
that $ is ill-defined or divergent when r vanishes. Note that this incorrect behavior of $ for 
D- 7^ is not a problem for the existence of a Lyapunov functional, but comes rather from 
the limited validity of the hypothesis used for its approximate construction. Nevertheless, as 
soon as the strong gradients disappear $ relaxes monotonously to the value $ = —0.79997, 
corresponding to the final state, a TW of wavenumber k = —0.0123. 

As another test in the non-chaotic region, for = — 1 and b r = 0.5 (c\ = — 1, c<i = —2) 
we use as initial condition an Eckhaus- unstable TW (k = 0.54 > kE = 0.48) slightly 
perturbed by noise. The system evolves to an Eckhaus-stable TW (k = 0.31) by decreasing 
its winding number (initially v = 44 and finally v = 26). Fig. ^ shows the evolution of $ from 
its initial value $(0) = —1.485 the final one $ = —1.77. Although there is a monotonously 
decreasing baseline, sharp peaks are observed corresponding to the vanishing of r associated 
with the changes in v. When v finally stops changing, so that A is close enough to the final 
TW, $ relaxes monotonously as in Fig. [|. 

It was explained in Sect. |TXT| that there are parameter ranges in which $ is smaller near 
the boundaries for existence of TW, that is near k = ±^Ja/b r , than for the homogeneous 
TW: k = 0. This happens for example for Di = 1, b r = 1.25 (ci = 1, c 2 = —0.8). The 
corresponding function is shown in Fig. || If this prediction is true, and if $ is a correct 
Lyapunov functional, evolution starting with one of these extreme and Eckhaus- unstable 
TW would not lead to any final TW, since this would increase the value of the Lyapunov 
functional. This would imply the existence for this value of the parameters of an attractor 
different from the TWs perhaps related to the Spatio- Temporal Intermittency phenomenon. 
We use as initial condition at the parameter values of Fig. an unstable TW of wavenumber 
k = 0.64 ($ w —0.81), slightly perturbed by noise. From Fig. |8|, the system should evolve 
to a state with a value of $ value even lower than that. What really happens can be seen 
in Fig. |9|. The system changes its winding number from the initial value v = 52, a process 
during which $ widely fluctuates and is not a correct Lyapunov functional, and ends-up in a 
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state of v = 5, with a value of $ larger than the initial one. After this the system relaxes to 
the associated stable TW of k = 2tiu/ L = 0.061 < k E = 0.23. As clearly stated by Graham 
and coworkers, the expressions for the potential are only valid for small gradients. Since k 
is a phase gradient, results such as Fig. || can only be trusted for k small enough. 

Finally, we show the behavior of $ in the Spatio- Temporal Intermittency regime. Since 
v is constantly changing in this regime it is clear that (|3.7| ) will not be a good Lyapunov 
functional and this simulation is included only for completeness. We take D{ = and 
b r = 0.5 (ci = 0, c 2 = —2) and choose as initial condition a TW with k = 0.44 > k E = 0.30 
($ = —1.89814), with a small amount of noise added. The TW decreases its winding number 
and the system reaches soon the disordered regime called Spatio- Temporal Intermittency. 
Fig. [l(| shows that the time evolution of $ is plagued with divergences, reflecting the fact 
that v is constantly changing (see inset). It is interesting to observe however that during the 
initial escape from the unstable TW $ shows a decreasing tendency, and that its average 
value in the chaotic regime, excluding the divergences, seems smaller than the initial one. 



V. NUMERICAL STUDIES OF THE LYAPUNOV FUNCTIONAL IN THE 

PHASE TURBULENCE REGIME 



The Phase Turbulence regime is characterized by the absence of phase singularities (thus 
v is constant). This distinguishes it as the only chaotic regime for which $ would be well- 
defined. Graham and co-workers [ 14 ,30f] derived especially for this region an expression 



proposed as Lyapunov functional in the small gradient approximation (|3.12p. 



We recall that the calculations in []14| , |30| predict that the phase turbulent attractor lies 



on a potential plateau, consisting of all the complex functions satisfying fl3.14p , in which all 
the unstable TW cycles are also embedded. The value of the potential on such plateau can 
be easily calculated by substituting in (|3.12|) an arbitrary TW, and the result is 



a 2 



*„ = -j- • (5.1) 



'T 



We note that this value does not depend on Di nor D r and then it is independent of Ci, the 
vertical position in the diagram of Fig. [3], within the phase turbulence region. 

In this section we take also a — 1. We perform different simulations for Di = 1.75 
and b r = 1.25 (ci = 1.75, c 2 = —0.8). In the first one, we start the evolution with the 
homogeneous oscillation solution (TW of k = 0). This solution is linearly unstable, but 
since no perturbation is added, the system does not escape from it. The potential value 
predicted by (|5.1[) is $ p z = —0.8. This value is reproduced by the numerical simulation up 
to the sixth significant figure for all times (Fig. [11], solid line). This agreement, and the fact 
that the unstable TW is maintained, gives confidence in our numerical procedure. 

In a second simulation, a smooth perturbation (of the form fie tqx with q = 0.049 and 
/i = 0.09) is added to the unstable TW and the result used as initial condition. This choice 
of perturbation was taken to remain as much as possible within the range of validity of the 
small gradient hypothesis. After a transient the perturbation grows and the TW is replaced 
by the phase turbulence state (the winding number remains fixed to 0). The corresponding 
evolution of $ is shown in Fig. |TT| (long-dashed line). The value of the potential increases 
from $ p / to a higher value, and then irregularly oscillates around it. Both the departure and 
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the fluctuation are very small, of the order of 10~ 4 times the value of $. Simulations with 
higher precisions confirm that these small discrepancies from the theoretical predictions are 
not an artifact of our numerics, but should be attributed to the terms with higher gradients 
which are not included in ( 3.12|) . As a conclusion, the prediction that the phase turbulence 
dynamics, driven by non-relaxational terms, maintains constant $ in a value equal to the 
one for TW is confirmed within a great accuracy. 

It is interesting however to study how systematic are the small deviations from the 
theory. To this end we repeat the launching of the TW with a small perturbation for several 
values of Di = ci, for the same value of b r as before. The prediction is that $ should be 
independent of C\. The inset in Fig. |11| shows that the theoretical value $ p z = —0.8 is 
attained near the BF line, and that as c\ is increased away from the BF line there are very 
small but systematic discrepancies. The values shown for the potential are time averages of 
its instantaneous values, and the error bars denote the standard deviation of the fluctuations 
around the average. 

Again for C\ = 1.75, C2 = —0.8, we perform another simulation (Fig. [11], short-dashed 
line) consisting in starting the system in a random Gaussan noise configuration, of amplitude 
0.01, and letting it to evolve towards the phase turbulence attractor. As in other cases, there 
is a transient in which $ is ill-defined since the winding number is constantly changing. 
After this $ decreases. This decreasing is not monotonous but presents small fluctuations 
around a decreasing trend. The decreasing finally stops and $ remains oscillating around 
approximately the same value as obtained from the perturbed TW initial condition. The 
final state has v = —1, so that in fact the attractor reached is different from the one in 
the previous runs [y = 0) but the difference is the smallest possible and the difference in 



value of the associated potentials can not be distinguished within the fluctuations of Fig. 1 1 



These observations confirm the idea of a potential which decreases as the system advances 
towards an attractor, and remains constant there, but at variance with the cases in the 
non-chaotic region here the decreasing is not perfectly monotonous, and the final value is 
only approximately constant. 

Since the small discrepancies with the theory increase far from the BF line, and since it 
is known that condition (|3.14j ) can be obtained from an adiabatic-following of the modulus 
to the phase that losses accuracy far from the BF line, one is lead to consider the role of 
adiabatic following on the validity of $ as a potential. To this end we evaluated <3> along 
trajectories A(x, t) constructed with the phase obtained from solutions of (|1.1| ), but with 
modulus replaced by ( |3.14| ), so enforcing the adiabatic following of the modulus to the 
phase. No significant improvement was obtained with respect to the cases in which the 
adiabatic following was not enforced since that, in fact, adiabatic following was quite welll 
accomplished by the solution of ( |1 . 1|) . Then it is not the fact that the solutions of (|1 . 1| ) do 



not fulfill (|3.14j ) exactly, but the absence of higher gradient terms in both (|3.14j ) and (|3.12|) 
the responsible for the small failures in the behavior of $. 

Finally, it is interesting to show that the Lyapunov potential $ can be used as a diagnostic 
tool for detecting changes in behavior that would be difficult to monitor by observing the 
complete state of the system. For example the time at which the phase turbulence attractor 
is reached can be readily identified from the time-behavior of $ in Fig. [II]. More interestingly 
it can be used to detect the escape from metastable states. For example, Fig. shows $ for 
evolution from a Gaussian noise initial condition (e = 0.01). D { = 2 and b r = 1.25 (ci = 2, 
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C2 = —0.8). The system reaches first a long lived state with v = 2 not too different from 
the usual phase turbulent state of v = 2. After a long time however the system leaves this 
metastable state and approaches a more ordered state that can be described as phase 
turbulent fluctuations around quasiperiodic configurations related to those of p8| . More 
details about this state will be described elsewhere ||36|| . What is of interest here is that 



from Fig. [12] one can easily identify the changes between the different dynamical regimes. 
In particular the decrease in the fluctuations of $ near t ~ 1000 identifies the jump from 
the first to the second turbulence regimes. 



VI. CONCLUSIONS AND OUTLOOK 

The validity of the expressions for the Lyapunov functional of the CGLE found by Gra- 
ham and coworkers has been numerically tested. The most important limitation is that they 
were explicitly constructed in a approximation limited to small gradients of modulus and 
phase. This precludes its use for evolution on attractors such that zeros of r and thus phase 
singularities appear (defect turbulence, bi-chaos, spatio-temporal intermittency) . The same 
problem applies to transient states of evolution towards more regular attractors, if phase sin- 
gularities appear in this transient (for instance decay of an Eckhaus unstable TW, evolution 
from random states close to A = 0, etc.). A major step forward would be the calculation 
of the Lyapunov potential for small gradients of the real and imaginary components of A, 
which would be a well behaved expansion despite the presence of phase singularities. 

Apart from this, if changes in winding number are avoided, expressions ( |3.6|) , ( |3.7|) , and 
( |3 . 1 2| ) display the correct properties of a Lyapunov functional: minima on stable attractors, 
where non-relaxational dynamics maintains it in a constant value, and decreasing value 
during approach to the attractor. These properties are completely satisfied in the non- 
chaotic region of parameter space, even in complex situations such as TW competition, 
as long as large gradients do not appear. It is remarkable that, although the potential is 
constructed trough an expansion around the k = TW, its minima identify exactly the 
remaining TW, its stability, and the non-relaxational terms calculated by substracting the 
potential terms to ( |1.1| ) give exactly their frequencies. In the phase turbulence regime, 
however, there are small discrepancies with respect to the theoretical predictions: lack of 
monotonicity in the approach to the attractor, small fluctuations around the asymptotic 
value, and small discrepancy between the values of the potential of TWs and of turbulent 
configurations, that were predicted to be equal. All these deviations are very small but 
systematic, and grow as we go deeper in the phase turbulence regime. They can be fixed in 
principle by calculating more terms in the gradient expansion. 

In addition in order to clarify the conceptual status of non-relaxational and non-potential 
dynamical systems one can ask about the utility of having approximate expressions for 
the Lyapunov functional of the CGLE. Several applications have been already developped 
for the case in which (|1 . 1|) is perturbed with random noise. In particular the stationary 
probability distribution is directly related to $, and in addition barriers and escape times 
from metastable TW have been calculated |12]J30|] . In the absence of random noise, $ should 
be still useful in stating the nonlinear stability of the different attractors. In practice however 
there will be limitations in the validity of the predictions, since $ has been constructed in 
an expansion which is safe only near one particular attractor (the homogeneous TW). 
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Once known powerful statistical mechanics techniques (mean field, renormalization 
group, etc. ) can in principle be applied to it to obtain information on the static properties 
of the CGLE (the dynamical properties, as time-correlation functions, would depend also on 
the non-relaxational terms N, as in critical dynamics ||15|| ). Zero-temperature Monte Carlo 
methods can also be applied to sample the phase turbulent attractors, as an alternative to 
following the dynamical evolution on it. All those promising developments will have to face 
first with the complexity of Eqs. ( |3.6p , ( |3.7| ), and ( |3.12| ). Another use of Lyapunov poten- 
tials (the one most used in equilibrium thermodynamics) is the identification of attractors 
by minimization instead of by solving the dynamical equations. In the case of the TW 
attractors, solving the Euler-Lagrange equations for the minimization of $ is in fact more 
complex than solving directly the CGLE with a TW ansatz. But the limit cycle character of 
the attractors, and their specific form, is derived, not guessed as when substituting the TW 
ansatz. For the case of chaotic attractors (as in the phase turbulence regime) minimization of 
potentials can provide a step towards the construction of inertial manifolds. In this respect 
it should be useful considering the relationships between the Lyapunov potential of Graham 
and coworkers and other objects based on functional norms used also to characterize chaotic 



attractors f37|,|38 
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FIGURES 



FIG. 1. Regions of the parameter [c\ — C2]-space (a = 1) for the d = 1 CGLE displaying 
different kinds of regular and chaotic behavior. Two analytically obtained lines, the Benjamin- Feir 
line (B-F line) and the D_ line, are also shown. 

FIG. 2. Relaxation to the simple attractor for a < 0. The parameter values are 
a = — 0.01,ci = 1 and C2 = —0.8 . The initial condition is a TW of arbitrary wavenumber 
k = 0.295 and arbitrary amplitude Aq = 1.0 . 

FIG. 3. Time evolution of <I> on the Z)_ line. The parameter values are 
a = 1, c\ = —1 and C2 = —1. The initial condition is a a Gaussian noise of amplitude e = 0.01. 
The system evolved towards a TW attractor of wavenumber k = 0.0245. 

FIG. 4. Time evolution of in the non-chaotic region for c\ = — 1 and C2 = —0.8 . The initial 
condition is an Eckhaus stable TW of wavenumber k = 0.13 perturbed by random noise of small 
amplitude e = 0.09. 

FIG. 5. Same as Fig. [I| but for c\ = —1 and C2 = —2 . The initial condition for A consists of 
two Eckhaus stable TW of different wavenumbers {k\ = 0.4, k<i = 0.08) joined together. The inset 
shows the real part of this initial configuration. 

FIG. 6. Same as Fig. || but for c\ = 1 and C2 = —0.8. The initial condition is a random noise 
of amplitude e = 0.01. 

FIG. 7. Same as Fig. || but for c\ = — 1 and C2 = —2 . The initial condition is an Eck- 
haus-unstable TW (k = 0.54 > kE = 0.48) slightly perturbed by noise. 

FIG. 8. The function ^ = 3>[-Afc] as a function of k. The parameter values are 
a = 1, c\ = 1 and C2 = —0.8 . The values of kE are indicated by dashed lines. The diamont 
indicates the point $^=0.64 taken as initial condition for the simulation in Fig. |9| 

FIG. 9. Time evolution of $ for c\ = 1 and C2 = —0.8. The initial condition is an Eck- 
haus-unstable TW {k = 0.64 > kE = 0.48) slightly perturbed by noise. 

FIG. 10. Time evolution of $ in the STI region (ci = 0.0 and C2 = —2). The initial condition 
is an Eckhaus-unstable TW {k = 0.45 > kE = 0.30) slightly perturbed by noise. The winding 
number evolution is plotted in the inset. 
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FIG. 11. Time evolution of $ in the Phase Turbulence region {c\ = 1.75 and C2 = —0.8). 
Solid line: evolution of a unperturbed unstable traveling wave. Dotted line: evolution from noise. 
Dashed line: evolution from a slightly perturbed traveling wave. The inset shows final average 
values of $ as a function of the c\ parameter (c2 = —0.8). The error bars indicate the standard 
deviation of the fluctuations around the average value. 



FIG. 12. Same as Fig. [TT| but for ci 
noise with an amplitude e 



2 and C2 



-0.8 . The initial condition was random 



0.01, time step 0.005. In this case 2048 Fourier modes were taken into 
account. Note the transition occurring arround t ~ 1000 to a less fluctuating state. 
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